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Abstract 

The properties of rapidly quenched superfluid phase transitions are computed for two-dimensional Kosterlitz- 
Thouless (KT) systems. The decay in the vortex-pair density and the recovery of the superfluid density after a 
quench are found by solving the Fokker-Planck equation describing the vortex dynamics, in conjunction with the 
KT recursion relations. The vortex density is found to decay approximately as the inverse of the time from the 
quench, in agreement with computer simulations and with scaling theories. 
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The properties of quenched phase transitions 
are of interest because they may be relevant to 
the rapidly cooling early universe [1]. The cosmic- 
string phase transitions in that case have recently 
been linked with the vortex-loop superfluid transi- 
tion in liquid helium and high-T c superconductors 
[2]. In this paper we investigate a new technique 
for calculating the quenched superfluid transition 
in two dimensions, employing Kosterlitz-Thouless 
renormalization methods. The results are in agree- 
ment with previous studies that used either com- 
puter simulations [3] or scaling theories [4] . 

The stochastic dynamics of the vortex pairs of 
the Kosterlitz-Thouless theory are modeled by a 
Fokker-Planck equation [5] for the distribution 
function T(r,t), which is the density of vortex 
pairs of separation r : 
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where I = ln(r/a ), with a the vortex core ra- 
dius, and U is the interaction potential between a 
pair. The time here is in units of the diffusion time 
of the smallest pairs of separation a , t' — t/r , 
where t^ 1 = 2D /a 2 , with D the diffusion coeffi- 
cient. From the KT theory 
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k B T dl 

wherein = h cr s /m 2 fcBTisthedimensionlessareal 
superfluid density. K is renormalized by the pres- 
ence of the vortex pairs, and the KT recursion re- 
lation for K can be written in terms of T as 
dK 

~m 

To quench from Tkt to low temperature it is 
first necessary to compute K (I) and the vortex dis- 
tribution at Tkt- The equilibrium solution of Eq. 
(1) is just the Boltzmann distribution, and taking 
the normalized core energy of a pair to be 7r 2 _R'(0), 
iteration of the KT recursion relations gives the 
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Fig. 1. Distribution function T as a function of the pair 
separation for different times after the quench. 
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Fig. 2. Vortex density and supcrfiuid fraction as a function 
of time. 

t' = vortex distribution shown in Fig. 1. The it- 
erations are taken to a maximum pair separation 
of I = 8 (r/a = 3000). 

We then (instantaneously) change the temper- 
ature to 0.1 T K t (K(0) -> lO-K^O)) and compute 
the new renormalized values of K{1) from Eq. (3), 
using a Rungc-Kutta method. Eq. (1) is then solved 
for the changed T after a time step At' with a two- 
step Lax-Wendroff algorithm. The changed vortex 
distribution is used to recompute K(l) from Eq. 
(3) at the increased time, and the entire sequence 
is continuously repeated to step out in time. In 
quenching to a low temperature the second term 
in the parentheses of the right side of Eq. (1) is the 
dominant term, and as a first approximation we 



have neglected the first term. 

The results for the pair distribution function (in 
units a~ 4 ) are shown in Fig. 1. This is a classic 
phase-ordering system, and as expected the small- 
est pairs decay away first, while the largest pairs 
remain "frozen" until the longest times. By inte- 
grating the distribution function the vortex density 
can be computed, shown in Fig. 2 (in units a~ 2 ). 
At long times the density decays approximately as 
t' _1 , in agreement with previous simulations and 
scaling theories [3,4], although it appears to be ap- 
proaching this quite slowly (the exponent is -1.17 
at t' = 1 and -1.12 at t' — 250), probably because 
K and the supcrfiuid fraction (K(8)/K(0), dashed 
curve in Fig. 2) are still changing fairly rapidly with 
time. The 1/t' behavior can be traced to the non- 
linear term proportional to T 2 that results from 
the substitution of Eqs. (2) and (3) into Eq. (1). 

A key feature of the present results is the scal- 
ing of the vortex decay with the diffusion time r , 
something generally neglected in previous treat- 
ments [1,3,4]. Estimating for helium films D ~ 1 • 
10~ 4 cm 2 /sanda o ~ 1-10~ 8 cm gives r ~ 1-10~ 12 
s, in agreement with estimates of the same quantity 
in three dimensions [6] . The very fast decay times 
may well explain the inability of millisecond-scale 
experiments to detect the vortices [7] . 
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